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Abstract 

The method of separation of variables is shown to apply to both the classical and 
quantum Neumann model. In the classical case this nicely yields the linearization of 
the flow on the Jacobian of the spectral curve. In the quantum case the Schrodinger 
equation separates into one-dimensional equations belonging to the class of generalized 
Lame differential equations. 
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1 Introduction 



The classical Neumann model describes the motion of a particle constrained to lie on an (A^ — 1)- 
sphere in N dimensional space and submitted to harmonic forces Originally this model has 
been solved by Neumann |l|] for = 3 by separation of variables in the Hamilton- Jacobi 
equation, and subsequently K. Uhlenbeck has shown that the model is Liouville integrable by 
exhibiting (A^ — 1) conserved quantities in involution 0. The classical solution of J. Moser 0] 
was shown to be deeply related to the theory of hyperelliptic curves by Mumford @] emphasizing 
the algebraic integrability of the model. 

Recently, new methods for solving integrable models have been introduced H, ^ ||, 
based on underlying Lie algebraic structures. The relation between the group-theoretical and 
the algebraic-geometrical approaches has been explained in . 

However the old methods of resolution by separation of variables begin to attract new inter- 
est As a matter of fact the solution of the quantum Toda chain was given in Gutzwiller [12 
and Sklyanin [O by a method closely connected to this procedure. 



We show in this article that the general Neumann model (any A^) admits such a resolution 
both at the classical and quantum level. This is achieved by choosing appropriate coordinates 
on the sphere, generalizing Neumann's coordinates for A^ = 3 |^. Considering what we 



know about the Lax formulation of the model |]14|, this choice of coordinates is also strongly 



reminiscent of the one advocated by Sklyanin [|11|]. The simplifying feature of this model as 



compared to the general case, is that the separation of variables is achieved by a transformation 



involving only the space coordinates and not the momenta |TT 

We first recall how the separation of variables in the Hamilton- Jacobi equation works in the 
classical case. This is a straightforward generalization of the work of Neumann ^. One gets 
(A^ — 1) one-dimensional ordinary differential equations, which are readily solved. We will see 
that the resolution of the Hamilton-Jacobi equation in this situation of a Liouville integrable 
system provides a canonical transformation which allows to linearize at one stroke the flows of 
the complete family of hamiltonians in involution. The separation of variables naturally lead 
to Abel sums and provide a simple proof of the linearization of the flow on the Jacobian of the 
Riemann surface associated to this dynamical problem by the Lax spectral theory P, |, [15 



In the quantum case, the Schrodinger equation similarly separates itself into (A^ — 1) one- 
dimensional equations. The quasi-classical quantization conditions are readily obtained. The 
implications of the exact quantization conditions will be studied elsewhere. 



2 The Neumann model 

The Neumann model has two equivalent formulations, one in terms of constrained dynamical 
variables with Dirac brackets, and an unconstrained one with canonical Poisson brackets. These 



formulations remain equivalent at the quantum level |]T5| so we shall stick to the unconstrained 
one which is more convenient for our present purposes. 

We start from a 2A^-dimensional free phase space {xn,yn, n = 1- ■ ■ N} with canonical 
Poisson brackets: {xn,ym} = ^nm and introduce the "angular momentum" antisymmetric ma- 
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trix: Jki = XkVi — xii/k and the Hamiltonian: 

^ k^l ^ 

We shall assume in the following that: ai < 02 < ■ ■ ■ < a^. The hamiltonian equations are, 
with X = (xk) Y = {t/k) and the diagonal constant matrix A = (akSki)'- 

X = -JX Y = -JY - AX 

They automatically ensure that J2xl remains constant and lead to the non-linear Newton 
equations for the particle: 

ifc= —(^kXk — Xk E(-^'^ ~ Q-i^f) 
I 

Let us note that the above hamiltonian is invariant under the canonical transformations: 

{xk, Vk) — > {xk, Vk + 4>Xk) (2) 

where may be any function of '}Zk^\- Classically the constrained system can be viewed as 
the hamiltonian reduction of the above unconstrained one under this symmetry. 

The Liouville integrability of this system is a consequence of the existence of (A^ — 1) 
independent quantities in involution, first found by K. Uhlenbeck: 

Fk = xl + Y.^^ (3) 

i^k ^fc ~ 

Notice that H = 1/2 Y,k o-kFk and Y,k Fk = ^■ 

3 Neumann's coordinates on 

Generalizing Neumann's choice in [|I|, we shall introduce (A^ — 1) parameters ti, ■ ■ ■ ,tAf-i on 
the sphere. They are the roots of the equation: 

Notice that this equation is invariant by Xk — ^ Ax^ so that ti < t2 < ■ ■ ■ < ti^^i are indeed 
defined on the sphere. Conversely, by definition of the tj we have for x G S''-^"^^: 



u{t)= ''') \ =^ xi 



Ukit - flfc) Ui^kiak - ai) 

Considering the graph of u{t) it is easy to see that: 

ai < ti < a2 < ^2 < ct3 < • ■ ■ < tN-i < (In 

and we have a bijection of this domain V of the tj's on the "quadrant" > V/c of the sphere. 
In general the sphere appears as a 2^-fold covering of the domain V, ramified on the edges of 
V. For instance Xi = y/ti — ai ip(t2, ■ ■ ■ , t^-i) so that xi changes its sign when we turn around 
the ramification point ti = ai in the complex plane. By similar analytic continuations we can 
cover the whole sphere. 
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Proposition. The (tj) 's define an orthogonal system of coordhiates on the sphere. 
A simple proof is obtained by considering, for each root tj the vector: 

^ f Xi Xn \ 

Vj — , ■ ■ • , 

\tj — Oi tj — I 

We easily check that x ■ Vj = 0, Vj ■ Vji = for j ^ j' and therefore {vj) is an orthogonal 
basis to the tangent space to S^^^^^ at x. On the other hand, taking the derivatives of the 
equations u{tj) = with respect to ti, - ■ ■ ,tAr_i, where the x^ are considered as functions of 
these variables, we obtain: 

Therefore: 

As a byproduct, since = —u'{tj) we get the metrical tensor: 

_ dx dx _ _1 Un^jitj - tn) 

4 The Hamilton— Jacobi equation 

It is convenient to start from the Lagrangian formalism since this is well suited to the change 
of coordinates on the sphere. 

1-2 1 

-U^-Y, 9jj'tjtj' - U 

33' 

Considering the polynomial: 

- au) u{t) ^J2 4 Hit - ai) = l[{t - tj) 

k k Ij^k j 

we see that: 

^ = ^ S ^kxl ^^(^ ak-J2 tj 
The conjugate momentum to tj is: Pj = dC/dij = Qjjij leading to the hamiltonian: 

H^Y.P3i3-C^-\Y.9"p]+U 

3 ^ 3 

where g^^ = {g~^)jj = ^/9jj- The Hamilton-Jacobi equation is a first-order non-linear partial 
differential equation obtained by substituting in H: pj dS/dtj. The action S is function of 
the space coordinates tj's. For a fixed energy E the Hamilton-Jacobi equation reads: 
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The method of separation of variables consists in looking for a so-called complete solution 
of eq. (^), i.e. depending on {N — 1) arbitrary constants, of the form: 

S(ti, ■ ■ ■ , tAT-i) = 5*1(^1) + ■ ■ ■ + S'Ar_i(tAr_i). 

For this purpose we need to consider the following Vandermonde determinant: 



1 


1 


tl ■ 




N-2 
1 


.N-2 



We shall also need similar determinants Dj obtained by removing the j column and the last 
row in D. One has the useful identities: 

D 



(6) 



Moreover considering the determinant obtained by replacing the last row by tj and expand- 
ing it over this row, we find: 



k = 0,---,N -3 

D k = N -2 
(Ei tj)D k = N-l 



(7) 



Finally denoting A{t) = nfc(^ — O'k) we have in view of 



9" = -^-If-'-' ^ Hti. 



and we can rewrite equation (H) as: 

-4 E (-1)^-^-^ D, A(t,) + D{- E + E - m = 

j \ j ) j k 



J \ J / 3 

Using identities (0) one can separate the variables, getting: 



4 A(tj- 



2 N-l 



+ E cfct) = 



fc=0 



where: 



cn-1 = 1, cn-2 = '2E — E (^k, and the other Ck arbitrary. 

k 



Let us remark that all S/s verify the same equation which we can write as: 



4A(t) 



dt 



2 AT-l 



(8) 



n=l 



where the 6„'s are (A^ — 1) independent constants. The energy is obtained in terms of the 6„'s 
as: E = 1/2 (J2k «fc - En K)- 
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5 The general linear flow on the Liouville torus. 



We have obtained the complete solution of the Hamilton- Jacobi equation: 

S{ti, ■ ■ ■ , t]sf-i; bi, - ■ ■ , ^AT-i) = E ^ji^j'^ bi, - ■ ■ , 6iv_i). 

j 

One can interpret this action as a generating function for a canonical transformation from the 
original variables {tj,Pj) to new variables ipn) linearizing the equations of motion. As usual, 
one has: 

_dS I _ dS 

This transformation is at the level of the symplectic structure and does not refer to any specific 
hamiltonian. Since all the in eq. play the same role it is natural to consider the Hamilton- 
Jacobi equations with respect to any one of the hamiltonians in involution F^. We show in this 
section that they all lead to the same separated equations for S. 

In order to do that, we first express Jki in terms of the coordinates tj and their conjugate 
momenta pj. Consider the canonical 1-form a = J^kVkdxk which reduces on the sphere to 
a = J2jPjdtj. So we have in view of eq. (H) 

Pj = 2^ yk-^ whence y ■ vj = 2pj . 
k ^^J 

This determines y modulo a vector proportional to x which does not affect the value of Jki- We 
find y = 1/2 J2j 9'''^Pj^j ^^d- 

Jki = -\ E (^fe ~ ^i) ^i^i 3"Pj 
With the help of this formula, Fk may be expressed as: 



It is convenient |]14| to introduce the generating function for the F^. 

W(A)-Et^ = ##^ (9) 

k ^-cik llfc(A-afc) 

for appropriate 6„'s and we have used J^^k = I- By a simple calculation we find: 

ni\)=ui\) (|i-E|^j (10) 

where u{X) has been introduced in section 3. Let us consider the Hamilton- Jacobi equation 
associated to the hamiltonian ?i(A). It reads: 
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Proposition. For any k = 0, ■ ■ ■ , N — 2 we have: 

This is proven similarly to eq. (|^) by considering the determinant obtained by replacing the 
last row of D by tj/{tj — A). Using now the identity (|12D for A; = one can rewrite eq. (0) in 
the form: 



Dt,-X 



4A(t,) (^^^ + tf-^ (t, - A) + n(A - 







This leads to the separated equations: 

4 A(t,) (^l^y + n(A - 6„) + (t, - A) (tf-' + j: Q(A)t;.) = 

We look for a solution of the form: S = J2 Sj(tj) independent of A. This is indeed possible and 
is achieved by setting: 

1=0 ^ ~ A 

and we recover exactly eq. (^. In addition we have obtained an interpretation of the b^s in 
eq. (§): they are the roots of 7i(A). So we have the expression of K. Uhlenbeck's conserved 
quantities Fk in terms of the integration constants 6„ of the Hamilton-Jacobi equation: 

P I\n{ak - bn) 

i'k 



Yli^ki^k - ai) 

We have shown that the 6„'s are equivalent to the F^'s as action variables. Since the 
symplectic form is: 

oj = da = dpj A dtj = ^ dbn A dipn 
j " 

we see that the 6„'s are in involution. This is still another proof of the involution property of 
the -Ffc's. Moreover their conjugate variables ipnS have a linear time evolution under Ti.{X) for 
any A. We have proven: 

Proposition. The complete solution of the family of Hamilton-Jacobi equations 7i(A) is 
J2j S(tj) for S solution of (j^ and ?i(A) is expressed in terms of the 6„'s only. 



6 Discussion of the classical solution. 

As an application let us solve the equations of motion for the coordinates tk for any hamiltonian 
7i(A). These equations are readily obtained: 

i, = {niX),t,} = -2u{X)^ with 



Introducing the polynomial: 



TV N-1 

Pit) = llit-ak) llit-bn 

k=l n=l 

and denoting by r the time variable we get: 

dtj 4:u{\)dT 



-P{tj) I{n^j{tj-tn){tj- X) 



(13) 



Now it is natural to introduce the hyperelliptic Riemann surface of genus (A^ — 1) defined 
in the complexified plane by the equation: 

+ p{t) = (14) 



As a matter of fact this curve also appears as the spectral curve, in the Lax pair approach . 
The independent abelian differentials of first class on this surface are the t^dt/s for k = 
0, ■ ■ ■ , N — 2. Recalling the identity ([12|) we get, for precisely these values of k: 

E ^ = 4drn(A) E (-1)"'^'^' = n n , (15) 
j Sj j tj- X D Hi (A - ai) 

Let us interpret geometrically this remarkable relation: to the point x on the sphere are 
associated (A^ — 1) unordered points (tj, sj) on the Riemann surface of equation (p!4D , the t/s 



being the roots of u{t) =0 and sj = y~P{tj). Such an unordered set defines a divisor on 
the surface, and since (A^ — 1) is the genus of this surface the path integrals of the Abel sums 
appearing in the left-hand side of equation ([T5|) provide an analytic bijection into the Jacobian 
torus of the Riemann curve. The right-hand side then shows that this point on the Jacobian 
evolves linearly with time. Of course the physical quantities must be real, so the divisor moves 
on a connected component of the real slice of the Jacobian, see |^ , which can be identified with 
the Liouville torus. 

Let us illustrate this geometry for Neumann's case of A^ = 3. The Riemann surface is of 
genus 2, and may be pictured as: 
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There are two abelian differentials: dt/s and tdt/s leading to two path integrals (Pq is any 
origin) : 

dt f(i2,s2) dt 
^1 = - + / - 

JPo S JPo S 

/•(ti.'^i) tdt /•(*2,«2) tdt 

^2 = + / 

JPo S JPo S 

Obviously Q = {Qi,Q2) has 4 periods corresponding to loops around a homology basis of the 
Riemann surface. The application {{ti, Si), {t2, S2)) — > ^ maps the Jacobian of the surface 
(analytic space of divisors modulo equivalence) to a 2-dimensional complex torus. 
Finally we notice that the action itself is given by a similar integral: 

2^ Jpo s 

Remark that the terms t'^ for = 0, ■ ■ ■ , A^ — 2 in Ylni'^ — bn) = + t'^ lead to abelian 

integrals of first kind, while the term t^~^ leads to an integral of second kind, having a double 
pole at 00. So the action may be seen as a multivalued meromorphic function on the Jacobian. 
More precisely, with similar notations, one can write: 

^ N-l 

•S" = - ^ Cfc-i + 5*0 
^ k=i 

where the c^'s are linear functions of F^s with coefficients polynomials in the a^'s, since Ylnit^ 
K) =Y{k{t- o-k) J2k Fk/{t - au) and 5*0 is: 

1 fPj t^~^ dt 

As a function of the divisor Pi + ■ ■ ■ + Pn-i, Sq has a simple pole on a variety of codimension 
1 obtained when one of the Pj's goes to 00, which is well known to be the divisor of a theta 
function on the Jacobian torus |Q. By restriction to real variables, S becomes a real multivalued 
analytic function on the {N — 1) -dimensional real Liouville torus. 



7 The quantum case 

The quantization of the Neumann system will be performed in the unconstrained formalism. 
The constraint will appear at the level of the wave functions. We consider 2N self-adjoint 
operators ?/„, n = 1, ■ ■ ■ , A^} with canonical commutation relations: 

[xk,xi] = [yk,yi] = 0, [xk,yi] = ih6ki . 

We define Jki = XkVi — xiUk so that Jkk = 0. Notice that there is no ordering ambiguity for 
k ^ I. Finally the quantum analog of K. Uhlenbeck's operators: 

Fk=xi+j:-^ 

are also defined unambiguously and are self-adjoint. 
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Proposition. The quantum theory is formally integrable, i.e. 



[Fk, Fi 







Proof. A proof of this fact was given in |]T^ using the i?-matrix formalism. Since this is a 
fundamental result we give here a more direct proof. One checks immediately that the Jki and 
Xk obey the algebra, similar to the classical Poisson algebra: 

[Jij, Jkl] = {^jkJli + ^ilJkj + ^jlJik + ^ikJjl) 



First one then shows that: 



E 



7^2 
•^kp 



We can assume k ^ I. Therefore in the commutator [J^p, J/g] we have three possibilities; p = q, 
q = k, and p = I. The case p = q produces the term: 



E 



7^2 



ai — On 







E 



1 



The case q = k gives: 

1 



\JkpJklJlp ~)~ JkpJlpJkl JklJlpJkp JlpJklJkp\ 



E 

p^i 



(ofc - ap)(afc - ai) 



\JkpJlpJlk ~l~ JkpJlkJlp ~l~ JlpJlkJkp ~l~ JlkJlpJkp\ 



Notice that this term vanishes for p = k so that one can assume p ^ k,l in the sum over p. 
Finally the case p = I similarly produces (with q ^ p): 



E 

p^l,k 



1 



{ai - ap){ak - ai) 



\JklJpkJlp ~l" JklJlpJpk ~l" JpkJlpJkl ~l" JlpJpkJt 



kl\ 



In this last expression, in the first and last terms one can bring Jki in the middle position. In 
doing so the commutators cancel. Then the three terms have the same J factors and add to 
zero. Similarly one shows that: 



X 



k 5 E 



q^l - 



X 



IT. 



'^kp 



p^k ^k Ctp 







so that finally [Fk, Fi] = 0. 

The canonical commutation relations are realized by setting: 



Vk 



-ih 



_d_ 

dxi 



We look for a simultaneous diagonalization of the F^s, i.e.: \E' = \E' . As in the classical 
case they are invariant under the symmetry (0) . This symmetry translates on the wave function 
by: 

vl> — > e*'^^'"^ \Ef . 
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Here is simply related to by = hjr if'{r) , and = J2k ^1- This shows that the symmetry 
is implemented by restricting oneself to wave functions independent of r. Consequently we will 
look for a wave function \E' depending on Neumann's variables tj defined on S^^~'^\ 

As a preparation to the general calculation, we consider first the Schrodinger equation 
associated to Neumann's hamiltonian: H = 1/2 J2k o>kFk- We find: 

^ k,i ^ fc 

It is easy to compute: 

J2 Jli = -h^ (2r^ + (4 - 2N) x.V - 2 {x.Vf ) 

k,l 

For wave functions independent of r we have x. V ^E' = and reduces to the orbital part 

of the Laplacian, i.e. the Laplacian on the sphere. So we can write without further ado the 
Schrodinger equation: 

(^^l\ + Uj^ = E^ (16) 
and express the Laplacian on the sphere using Neumann's coordinates as: 

Here g = Ylj gjj is the determinant of the metrical tensor. It involves Hj Unj^ji'tj — tn) = 
=tnm>n(^m — tn)^ = ±-D^ SO neglcctiug numerical factors that cancel themselves between ^ 



and 1 / y/g we can set ^/g = D/ JWn ^{tn) > 0. From this we get: 



and we remark that Dj/ JYln^j A(t„) is independent of tj and can be commuted with d/dtj to 
the left of A. 

Finally the equation (16) reduces to the following simple form: 



It is now easy to separate the variables, i.e. to look for a solution of the form: 

^(ti,---,t7v-i) = ^i(ti)^2(t2)---^jv-i(tjv-i) 
and we find that the \E'j's all obey the same ordinary differential equation: 



d /-—— d 



N-2 



Ah' V^A(t) J ^A{t) J - t^-' + (^ a, -2E)+Y: Ck t^-'-' ^{t) = 

k k k=l ) 
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Writing t^-i - (E. a, - 2E) - YZ=i t^''^"^ 

Y{n{t - hn) ^ 

Y\k{t -ak) ^ t-ak 
we get the separated Schrodinger equation in the form: 

1 



Unit - hn) and: 

fk 



(P 1 

— + -y 



d_ ^ fk 



t — ttk dt Ah? t — ak 



^{t) = 



(18) 



This equation has been studied in the htterature and is known as the generahzed Lame equa- 
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8 The general Schrodinger equation. 

Since all the F^'s commute and are self-adjoint, there exists a complete set of common eigen- 
vectors. We therefore expect that similarly to the classical case the same separation of variables 
occur for the whole family of hamiltonians ?^(A). We shall show that it is indeed the case. As 
a byproduct, this is an alternative proof that the FkS commute, and shows that the fkS in 
equation (|18D are their eigenvalues. 

For wave functions independent of r we have by (^), denoting Vj = Xk/{tj — dk)'- 



d 

dxk 



E 9" 



d_ 



ih 



Jki = ^-^Y.i(^k-ai)v^^v\g" 



Let us compute the generic hamiltonian: 

Fk 

k ^ ^ '^k 



d_ 



u 



(A) + iE 



■J I 



kl 



kl 



(A - afc)(A - ai) 



Proposition. We have (compare with the classical formula 



n{\) = u{\) 



^^E 



1 d ( 



V99 



X-tiy/gdti \ dti J 

Proof. We only sketch the main steps of the calculation. We first have: 



7^2 



y — 

^ (A - afc)(A - ai] 



4 ^ 

k,l,i,j 



{ak - CLi) 



(A - afc)(A - ai 



The summation over k can be performed with the help of the following identities: 

4 Qii^ij 



E '^^-^^ 



] ^ 



E 



k h 

X- ak 



4 ti Cjii^ij 



X-t 



u{X) 



(t,-A)(t,-A) 



97^ 
dti 
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leading to the result: 



nix) = uix)--Y: 



(A - ai)u{X) 4:gii{ti - aif5.. 



t-:A{U - X){tj - X) {U-X){ai-X)\ 



n33 



r — 



To perform the summation over I in addition to the previous identities we need the following 
ones: 



, dv] 
]■ — - 

1 , dvl 
> V — - 

^X-ai 'dU 

The following result then emerges: 



2^A.. + 2 ^" (1-6- 



ti tj 



ti tj 



ld_ 
2dt, 



u{X) gu 



(t,-A)2 X-U 



niX)=u{X) 



1 



1 <' 



d_ 

dt,. 



dti tj ti 



d_ 

dti 



The result now follows by noticing that: 

^logv^ _ 1 / iidgu 
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dti 2 

The corresponding Schrodinger equation is: 



-E 



dti tj ti 



n{X) ^=(j: 

\ k 



fk 



and defining the bnS such that Uni^ ~ bn) / ^ (A) = Y,k fk/{X — ak) it can be rewritten as: 



h D{X-t,) 



d 



d 



Ah' v/A(t,) — v/A(t,) — - tf-^ (t, - A) - n(A - v[; = 



As in the classical case this leads to the separated equations: 

\Ah' ^ I- - - bn) - {t, - A) ftf -2 + y Ck{X)t) ) \ v^,(t,) = . 



dti 



We obtain equations independent of A by choosing: 



Ar-3 

E 

fc=0 



n„(t - hn) - n„(A - 



i^-'+Ec.(A)t^ 



fc=0 



t-X 



and we get back eq. (0) with the /^'s identified with the eigenvalues of the -F^'s. 
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9 The semi— classical quantization conditions 



The solution of the quantum Neumann model has been reduced to the study of the Lame 



equation (18). We shall not dwell further on this question in this paper and content ourselves 
with a discussion of the semi-classical quantization. The semi-classical wave function is of the 
form: 

*/i-/iv(''^i'---'^iV-i) = exp (^^S{ti, - ■ ■ ,tN-i] 
Notice that the semi-classical quantization conditions: 



Pk dqk = 27rhn 



are canonical invariants and mean that \& is univalued on the Liouville torus since § pdq is the 
variation of S when one describes a non-trivial loop around the Liouville torus. 

In order to understand further these conditions, it is necessary to describe the corresponding 
cycles on the curve ([T^). Since P{t) = Y[k{'t~'^k) Y[n(J^'~^n) one has to worry about the relative 
disposition of the a^s and 6„'s which depends on the signs of the f^s. This introduces several 
cases as first noted by Neumann [jl|. In order to simplify the discussion we shall assume for 
example that all f^s are positive, implying: 

ai < bi < a2 < b2 < ■ ■ ■ < bN^i < o-n 

as may be seen immediately by considering the graph of J2kfk/{t — cLk)- Moreover we know 
that: 

tti < ti < a2 < t2 < < ■ ■ ■ < tj\f_i < Oat 

so that classically each tj is constrained in an interval [aj, bj] or [bj, ctj+i] in view of equation (H) 
i.e. each P(tj) must be < 0. Since P{t) > for t > ajy one sees that tj G [bj, flj+i]. 

Let us forget momentarily ti, ■ ■ ■ ,tN-i and concentrate on the motion of tj. Recalling 
equation (0) we see that dtj/dr vanishes when tj bj and tj aj+i and is of a fixed sign 
in this interval. So tj oscillates between its bounds. Now when tj = a^+i in view of section 
3, Xj+i = but Xj+i 7^ 0. So the point on the sphere S^^~^^ crosses the hmit Xj+i = 
of its quadrant when tj a^+i and continues its motion on the symmetric quadrant with 
respect to the hyperplane Xj+i = 0. In the space of the t variables this may be seen as an 
analytic continuation around the branch point tj = aj+i. To sum up our discussion one loop 
on the sphere corresponds to two oscillations on the interval [bj,aj^i]. Writing that the wave 
function resumes its initial value at the end of the loop we get the semi-classical quantization 
conditions: 

fOj + l 1 



— ^ dt = 2'iThni , n,- G Z 

2^ Uk{t-ak) ' 

These (A^ — 1) quantization conditions determine the 6„'s or equivalently the f^s. Of course 
the range of the /^'s should be restricted so that the previously mentioned constraints on the /^'s 
be obeyed (including in particular the positivity of the energy for the case of original Neumann's 
hamiltonian) . Finally let us notice that the integrals appearing in these conditions express 
themselves in terms of some of the periods of abelian integrals on the Riemann surface (|1^. 
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10 Conclusion. 



In this article we have shown that the method of separation of variables which was known to 
apply to the classical Neumann model applies equally well to the quantum Neumann model. 
This model is presumably one of the simplest non trivial integrable model which can be solved 
in this way. This is because the necessary change of variables to separate the equations is only 
a change of coordinates on the sphere. Moreover we have treated on the same footing the whole 
family of integrals of motion in involution discovered by K. Uhlenbeck in the classical case as 
well as in the quantum case. This method of resolution introduces a Riemann surface and the 
associated Jacobian torus in a particularly simple and elegant way. 

In the classical case this leads to the linearization of the flow on the Jacobian, while in the 
quantum case the quantization conditions involve the real periods of the underlying Riemann 
surface. We will promote our semi-classical analysis to an exact one in a coming work. 

Recently Sklyanin has reconsidered this method of separation of variables in various in- 
tegrable systems [|rT|. It is interesting to remark that in such system as the Toda chain the 
separated equations appear to be finite difference equations (Baxter equation) while in our case 
we end up with ordinary differential equations. We expect this to be related to the fact that 



in the Lax pair formulation the Neumann model involves a linear Poisson structure [14| while 



in the Toda chain the Poisson structure is quadratic ||T3 
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